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ABSTRACT 

In recent years, gravitational lensing has been used as a means to detect substructure in galaxy- 
sized halos, via anomalous flux ratios in quadruply-imaged lenses. In addition to causing anomalous 
flux ratios, substructure may also perturb the positions of lensed images at observable levels. In this 
paper, we numerically investigate the scale of such astrometric perturbations using realistic models of 
substructure distributions. Substructure distributions that project clumps near the Einstein radius of 
the lens result in perturbations that are the least degenerate with the best-fit smooth macromodel, 
with residuals at the milliarcsecond scale. Degeneracies between the center of the lens potential 
and astrometric perturbations suggest that milliarcsecond constraints on the center of the lensing 
potential boost the observed astrometric perturbations by an order of magnitude compared to leaving 
the center of the lens as a free parameter. In addition, we discuss methods of substructure detection 
via astrometric perturbations that avoid full lens modeling in favor of local image observables and 
also discuss modeling of systems with luminous satellites to constrain the masses of those satellites. 
Subject headings: cosmology: theory - dark matter - gravitational lensing 



1. INTRODUCTION 

The cold dark matter (CDM) paradigm predicts that 
the dark matter (DM) halos hosting most galaxies con- 
tain a large number of low-mass, compact subhalos 
within their virialized regions. These subclumps, collec- 
tively referred to as substructure, may or may not con- 
tain luminous stellar tracers. Since the observed number 
of dwarf galaxy satellites in the Local Group falls short of 
the expected n umber of subhalos by more than an orde r 
of magnitude (iKl vpin et aI1ll999bi : iMoore et al1ll999f ). 
then either the CDM model is incorrect or dwarf galax- 
ies are bia sed tracers of DM in gala xy-scale h alos (e.g. 



Spergel fc Stcinhardt 2000: Hanncst ad fc Sche rrer 2000; 
Hu et a l.' '2000: Bifllock et al.ji200a iBenson et al.ii2002i : 
Somerv iUc 2002; Stochif [200a iNagai fc Kravtsovl (20051 ). 



If the CDM paradigm is correct, then the paucity of op- 
tical counterparts of substructures leaves us with few av- 
enues for detecting and investigating DM substructure. 

One possibility for studying dark substructure 
around local galaxies is through its perturbative ef- 
fect o n kinematica ll y cold systems like tidal streams 
Ce.g.. iMaver et all 120021) or galactic disks (e.g., 
iToth fc Ostrikeilll992t iBenson et al.ll2004D . These tests 
probe only the closest and most massive substructures, 
and are only applicable in very nearby galaxies. As a 
result, this method provides little information to address 
fundamental questions such as the amplitude of the mat- 
ter power spectrum on subgalactic scales or the detailed 
properties of dark matter itself. 
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Another promising method of study is through the 
gravitational lensing effects of substructure, referred to 
as substructure lensing. Much of the previous work 
in this field has focused on the effec ts of substructure 
on image magnifications and fluxes (iMao fc Schneideij 
iggi; 'Metcaff fc Madaul [200H iDa lal fc Kochanekl l200a 



Kochanck fc Dalai 2004 jRozo et al. 2006) or on time de- 
lays ( Morgan et al.ll200a l. However, substructure can 
perturb the deflection angle, a, of lensed images as well. 
Image positions in lensed systems may not be subject to 
the same foregrounds - such as dust absorption - as flux 
anomalies. Further, since the astrometric perturbation, 
5a, is a steeper function of subhalo mass than flux per- 
turbations, it may provide a qualitatively distinct probe 
of substructure properties. 

A study of astrometric perturbations is particularly 
timely given the possibility of submilliarcsecond res- 
olution obs ervations of strong lenses. For example, 
iBiggs et al.l ((2004) have presented observations of a 
four image jet source system, CLASS B0128-I-437, in 
which milliarcseond astrometric perturbations may have 
been detected. Previous studies of astrometric per- 
turbations by substructures have come to conflicting 
conclusions about their overall size and probability. 
iMetcalf fc Madaul (|20f)lh use lensing simulations of ran- 
dom realizations of substructure in regions near images 
to suggest that in order to change image positions by 
a few tens of milliarcseconds (mas), there must be sub- 
clumps with masses > lO^M© in Milky Way sized ha- 
los that are well aligned with the images they perturb. 
This alignment is likely to be rare in CDM, although the 
probability would in crease in syste ms where the source 
is elongated in a jet ()MetcaHll2002[ ). On the other hand, 
Chiba (200^ tests the size of astrometric perturbations 
in B1422-I-231 with a model of CDM subhalos as point 
masses and finds deflections of 10 to 20 mas, using sub- 
clumps with masses greater than ^ 2 x lO^/i^^M©. Ad- 
ditional studies of astrometric perturbations by single 
perturbers have focused on the detecti on of substructure 
via distortion of a finite source (jlnoue fc Chibal l2003l . 
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l2005allH ). In addition, iPen fc Mad (|2005f) studied the 
effect of multiple lens planes on the rotation of lensed 
images. 

In this paper, we estimate the amplitude of astrometric 
perturbations produced by substructures using realistic 
substructure models and test the feasibility of observ- 
ing such perturbations by comparison with the image 
positions given by a be st-fit lens macromodel. While 
iMetcalf k, Madaul(|2001| ) suggest that single clumps need 
to be well aligned with the images in order to produce 
observable perturbations, measurable effects may be pos- 
sible in more general scenarios, given the collective ef- 
fects of entire substructure distributions. In addition, 
the kinds of perturbations in these scenarios may be dif- 
ferent. For example, while the scenario of a nearby sub- 
halo that perturbs only a single image ensures that de- 
generacies with best-fit smooth macromodels are small, 
large numbers of more distant substructure may affect 
multiple images and exhibit large degeneracies with the 
macromodel. In addition to using detailed lens model- 
ing, we discuss model-independent methods for identi- 
fying substructures, proposing a method using systems 
with multiply imaged jets. 

Substructure lensing may also be used for comparisons 
between luminous satellites and dark subhalos. For ex- 
ample, we can address one of the open questions regard- 
ing CDM substructure: why are some subhalos dark and 
some luminous? There have been many proposed mecha- 
nisms to explain which subhalos have stars and which do 
not. One possibility is that the efficiency of star forma- 
tion diminishes with decreasing halo mass, perhap s due 
to photoionization squelching (Bullo ck et al. I [20011 ). An- 
other suggestion has been that luminous galaxies form 
only in the highest mass halos {M > IO^Mq), and over 
time tidal stripping reduces subhalo masses to the low 
values (M ~ 1 0^ Alp,) inferred for the smallest local 
dwarfs (jKravtsov et al.ll2004bf l. Yet another possibility 
is that the low masses inferred for local dwarfs are in 
fact mistaken, and i nstea d are much larger, M >10^ Mq 
(jStoehr et al.ll200l |2003| ). iKazantzidis et all (|200i ) find 
that subhalos do not experience the significant mass re- 
distribution in their centers required to embed satellites 
in massive subhalos, and studies stellar kinematics in 
the d warfs also dispute the e xistence of massive subhalos 
(e.g.. [Wilkinson et al.|[200l . Here, we propose directly 
measuring the masses associated with luminous satellites 
projected near lensed images , the r esults of which would 
direcly test the IStoehr et al.l (|2002( ) hypothesis. 

The layout of the paper is as follows. We begin by de- 
scribing the host halo and substructure models we will 
be employing in our analysis in ^21 The statistical prop- 
erties of the astrometric perturbations in our models are 
presented in ^ and their degeneracy with macro-lens 
model parameters is discussed in ^ We discuss model- 
independent, astrometric signatures of substructure lens- 
ing in SjS] Finally, the effect of luminous substructures 
on astrometric perturbations and lens modeling is dis- 
cussed in 511 and we present a summary of our results 
and conclusions in section ^JT] 

2. NUMERICAL ESTIMATES 

In this section we present the host halo and substruc- 
ture models we used to analyze astrometric perturbations 
and describe our methods of creating artificial lens sys- 



tems with substructure, finding the image positions in 
such systems, and fitting for the best-fit smooth model. 

2.1. Halo Model 

We begin by specifying the macromodel parameters of 
the lens system - using typical lens parameters ~ and the 
background cosmology used to generate artificial lenses. 
In particular, we choose a lens redshift of zi — 0.5, a 
source redshift of Zs = 2.0, and a flat ACDM cosmology 
with flm = 0.3, Oa = 0.7, and h = 0.7. The halo of 
the lensing galaxy is modeled as a singular isothermal 
ellipsoid (SIE), the projected density profile of which is 
given by 

'^(O = ^, (1) 

where ^ is the projected radius, — x"^ + y'^/q^ and 
q = 0.9 is our fiducial value for the lens's axis ratio. 
The particular values for the ellipticity and shear do 
not appear to be significant to the results of this pa- 
per. The length scale b corresponds to the Einstein ra- 
dius of the lens for circularly symmetric profiles {q — 1). 
We take as our fiducial value b = 1.05", corresponding 
to the Einstein radius of a singular isothermal sphere 
of mass M — IO^^Mq in our chosen cosmology and 
source and lens redshifts. Finally, we assume there is 
an external tidal shear of 7 = 0.16 aligned with the 
major axis of the halo {6-y — 0). Our particular choice 
7 = 0.16 correspo nds to the best- f it valu e for the external 
shear obtained bv iBradac et al.l (|2002f ) for the lens sys- 
tem B1422-I-231, a system known to exhibit anomalous 
flux ratios. 

2.2. Substructure Models 

We employ four different subhalo catalogs to populate 
our fiducial halo with substructures. The catalogs de- 
scribe the three-dimensional position in the parent halo 
and the density profile of each subhalo. In order to 
use these catalogs in our halo model, the substructure 
positions and density profiles must be projected along 
an axis. The density profiles are projected by param- 
eterizing the spherical density profiles in three dimen- 
sions, then using an approximate formula for the profile 
in two dimensions. The following section describes the 
substructure models in detail. 

For two of these models, we use substructure catalogs 
from a DM cluster simulation with flm = 0.3, J^a = 
0.7, h = 0.7, and as = 0.9. The adaptive -refinement-tree 
(ART) code is used to run the simulation (jKravtsov et al.l 
flWTl: lKra^^tso^fl99^ using a 256'^ uniform root grid 
covering the computational box of 80h~^ Mpc and res- 
imulating a cluster in the box with more particl es and 
higher spatial resolution (T asitsiomi et al.ll2004D . For 
this higher resolution simulation, the effective mass reso- 
lution is rup — 3.95 x 10^h~^ Mq and the smallest cell size 
reached is 0.6h~^ comoving kiloparsecs. Detail s of th e 
simulations can be found in iNagai fc Kravtsovl ([2005). 
A variant of the Bound Dens ity Maxima halo finding 
algorithm (jKlvpin et al.lll999ah is used to identify sub- 
halos. Deta ils of the algorithm used to find subhalos can 
be found in iKravtsov et all ()2004af l. 

We draw two different subhalo catalogs from the clus- 
ter simulation, one at z=0.5, as the cluster is undergoing 
a major merger, and another at z=0. The mass of the 
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parent halo itself is scaled down to 10^'^ Mq and all sub- 
clumps are scaled accordingly^. Finally, clumps within 
the virial radius with mass M > IO^Mq {M/M^ir > 
10~^) are chosen. These clumps correspond to Np > 80. 
iNagai fc Kravtsovl ^2005) show that the mass function 
of the high resolution simulation converges with that of 
a lower resolution simulation for Np > 80. W e fit the 
clumps with a Moore profile (|Moore 



p{r) 



(2) 



(r/r,)i-5(l + (r/r,)i-5) 

truncated at a radius using a Levenberg-Marquardt 
minimization technique. The resulting overall mass frac- 
tion in substructure is 10% at z=0.5 and 13% at z=0. 
The left-hand panels of Figures [T] and [2] show the num- 
ber density profile and cumulative mass function for both 
redshift samples. 
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Fig. 1. — Left: The number density profile of subhalos from nu- 
merical models normalized to the number density within the three 
dimensional virial radius, n{r) / (n^ij.) , for z=0.5 (top) and z=0 
(bottom). The dashed line shows the dark matter profile normal- 
ized to the virial overdensity. Right: The number density profile 
of subhalos from semi-analytic models normalized to the number 
density within the virial radius, n(r)/(nvir), for a young halo (top) 
and an old halo (bottom), using all clumps with masses greater 
than M = 10* Mq. The dashed line shows a Moore profile with 
concentration of 6.4, normalized to the scale radius. 

Compared to galactic halo, the scaled cluster sim- 
ulation may have a greater number of substructures 
(jZentner et al.l20d5a[ l - although in our substructure cat- 
alogs, the substructure mass fraction does not seem sig- 
nificantly larger than expected at 10-20%. In addition, 
we may expect that the concentration of the subclumps 
in the cluster simulation may be smaller than in a com- 
parable galactic halo; the issue of subhalo concentrations 
is further discussed in Section !j4.4l 

The radial distribution of subhalos may be weakly 
dependent upon the mass scale of the parent halo 
(|Gao et all [20041 : iDiemand et"all I2004D . with a cluster- 
scale simulation with fewer subclumps in the center of the 

^ Although the two catalogs are drawn from different epochs, we 
use both for ad hoc substructure distributions in the lens at the 
lens redshift and, therefore, fix the masses and virial radii to the 
same values. 
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Fig. 2. — Left: The cumulative mass function of subhalos from 
numerical models enclosed within the virial radius, for subhalos 
at ^=0.5 (top) and 2=0 (bottom). Right: The cumulative mass 
function of subhalos from semi-analytic models enclosed within the 
virial radius for subhalos for a young halo (top) and an old halo 
(bottom). In all panels, the mass function from the three dimen- 
sional mass profile is a dotted line, the mass function derived from 
the parameterization of the subhalos using Eqn. |4] is a solid line 
and the normalized n ~ m~^ * mass function is a dashed line. 

halo than a galaxy-scale simulation. In addition, simula- 
tions may suffer from overmerging in the center of halos, 
reducing th e number of subclum ps in the inner po rtions 
of the halo (jMoore et ah 1996: Kl vpin et al.lll999af l. This 
effect may be particularly important in lensing studies, 
where substructure close to the projected center of the 
system produces the strongest perturbations. To get an 
independent estimate of how much substructure might 
survive in the inner parts of galaxy halos, we consider 
substructure distributions in two halos cons tructed us- 
ing the semi -analytic model described in Tavlor &: Babull 
(|200lL I2004f ): these models - like our simulation models 
- do not include a galaxy component, but contain much 
more substructure in the central regions of the cluster 
(|Tavlor fc Babull 120051) . thus likely providing an upper 
limit to the amount of substructure in CDM halos. 

The semi-analytic models assume spherical symmetry 
in the input orbits and the potential, and all clumps are 
parameterized by a Moore profile, trunc ated at rt and 
decreased in density by a fraction ft as in iHavashi et al.l 
(12003) (equation 8): 



p(r) 



ft 



(3) 



l + (r/rt)3 ((r/r,)i-5)(l + (r/r,)i-5)' 

'Kazantzidi s et ahl (|2004) suggests that the IHavashi et all 
(2003) relation may reflect lowered concentrations in the 
the subhalos simulation due to nonequilibrium initial 
conditions combined with numerical resolution rather 
than tidal shocks; the issue of subhalo concentrations 
is further discussed in Section ^14. 41 

The two semi-analytical substructure catalogs we 
choose correspond to a dynamically old and a dynam- 
ically young parent halo at z=0. The old halo has ac- 
creted 50% of the parent halo's mass by z=3.2 and 90% 
by z =0.67. The corresponding redshifts for the young 
halo are 0.54 and 0.04. Finally, the old halo has a sub- 
structure mass fraction of 10%, while the young halo has 
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a mass fraction of 21%. We scale the host halo mass - a 
galactic halo of LGxlO^^M© - to lO^^M© and scale ah 
the substructures accordingly, keeping only those sub- 
structures with rescaled masses larger than 10^ Mq. The 
right-hand panels of Figures [T] and [2] show the resulting 
number density profile and cumulative mass function for 
both samples. 

In comparison to the simulation profiles, the semi- 
analytic catalogs have many more clumps in the center of 
the halo, while the number density profiles of all the sub- 
structure models track the host halo dark matter profile 
at large radii. In addition, the mass functions of all the 
substructure models are similar, n m~^'*, as expected 
from the results of h igh resolution numerical simulations 
(|Ghigna et al.|[2000l ). In both the simulation substruc- 
ture models and the semi-analytic models, one of the 
distributions follows n ^ m~^ * to masses greater than 
1% of the halo mass, while the other distribution falls 
off and has fewer large-mass clumps. We shall see that 
the presence or absence of large mass clumps can have 
important consequences on the distribution of position 
perturbations. 

In addition to specifying the general properties of the 
host halo and its substructure population, lensing calcu- 
lations require we specify a line-of-sight projection axis. 
We have chosen to project all halos along the major axis 
of the host halo. This choice for line-of-sight projection 
is motivated by the fact that it results in the most com- 
pact, and therefore most effective, lenses. Consequently, 
one might expect most observed strong lens systems to 
be projected along this axis. Moreover, this projection 
leads to the largest substructure densities in the cen- 
tral part of the halo, which should maxim ize the impact 
substructures have on the lensed images (jZentner et al.l 
I2005H) . 

As expected, the projected substructure distributions 
obtained from the simulation catalogs have fewer clumps 
in the central regions of the halo than the projected semi- 
analytic models. For instance, within 2" of the center 
of the lens, the z = 0.5 simulation projection has no 
clumps and the z = simulation has 1 clump of mass 
3.8 X IO^Mq, or ~ 0.5% of the halo mass that is enclosed 
within the Einstein radius. On the other hand, both 
semi-analytic models have 7 clumps within this projected 
radius, with total masses of6.5xlO^Mo and T.lxlO^M© 
for the dynamically old and dynamically young catalogs 
respectively, or ~ 5 - 10% of the halo mass that is en- 
closed within the Einstein radius. For future reference, 
we note that the z = simulation model has the largest 
subclump of all the models, ~ 7% of the halo mass, and 
this subclump is projected far from the Einstein radius 
of the host halo. 

Rather than modeling the two dimensional mass profile 
of each individual substructure as a projected Moore or 
Hayashi profile, we use the approximation, 



where Kg — 2psrs for projected Moore profiles and 
Kg = 2ftPsrs for the Hayashi profile, and x = 2.AR/rs 
where R is the projected radius. Our choice of profile is 
numerically motivated in that the deflection angle asso- 
ciated with our chosen profile has a simple closed form 
expression, whereas the projected Moore profile does not. 



We show the mass functions of the substructure distri- 
butions using this approximation in Fig. [2] (solid lines) . 
Note our approximation does not appreciably bias the 
substructure mass functions. 

2.3. Generating and Modeling Artificial Lenses 

Our host halo and substructure models are used to 
generate artificial lenses as follows. First, the image po- 
sitions of a given point source are obtained by finding 
all roots of the lens equation using a Newton-Raphson 
root finder with a gridded set of initial guesses. Once 
the image positions of a source are obtained, we fit the 
resulting lens system with a lens model that parameter- 
izes the smooth components of the lens only. We then 
compare the true image position with that obtained in 
the absence of substructures and with that predicted by 
the best-fit model for each image in the lens. We repeat 
this experiment to statistically sample the source plane 
in order to obtain distributions for the astrometric per- 
turbations generated by the dark matter substructures. 

The results of such a statistical comparison obviously 
depend on the how the source plane is sampled. While 
naively one might expect uniform source plane sampling 
to be adequate, observed lens samples suffer from magni- 
fication bias (brighter systems are more likely to be ob- 
served than dimmer systems) and consequently magnifi- 
cation weighted sarn p ling is the most appropriate choice. 
iKeeton fc Zabludofj (|2004D show that for sources with 
a power law luminosity function, dN/dS cx S^"^ with 
V = 2, magnification weighting in the source plane is 
equivalent to uniform sampling of the image plane. The 
importance of this result rests on the fact that uniform 
image plane sampling is easy to implement and the for- 
tuitous fact that the observed quasar luminosity function 
in the largest lens surveys is roughly a power law with 
an index, j/ « 2. Consequently, we have opted for dis- 
tributing sources along the source plane in accordance 
with uniform image plane weighting in order to provide 
a closer match to observations. 

Best-fit models to artificial lenses are obtained by 
minimization with a downhill simplex algorithm. If con- 
vergence is not achieved within a prescribed number of 
steps, the original input parameters (the macromodel pa- 
rameters) are perturbed and the modeling is repeated 
until convergence is achieved. The minimization puts 
priors on parity agreement and excessively magnified im- 
ages during modeling. We eliminate extremely dim im- 
ages, fj, < 0.01, and systems with excessively magnified 
images, fi > 50. Finally, an appropriate best-fit model 
may not be found if the simplex travels to an area of pa- 
rameter space that produces a number of modeled images 
that is different from the observed number of images. 

The smooth models used in the minimization are com- 
prised of an SIE mass distribution along with an external 
shear component. The model parameters are the Ein- 
stein radius 6, the projected axis ratio q, the external 
shear 7, the orientation of the shear 9j, the orientation 
of the ellipticity 9q, the center of the potential ^origin and 
J/origin, f-nd the Coordinates of the source position Xgourco 
and ^source ■ 

Notice that lens systems with a single source position 
have nine parameters to be fit whereas there are only 
eight observables in quadruply imaged systems. To make 
the system overconstrained, we artificially hold the ori- 
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gin (i.e., the position of the SIE) fixed at {x,y) = (0,0), 
thereby reducing the number of free parameters to seven. 
ObservationaUy, this can be accomphshed by fixing the 
center of the mass distribution to the observed position 
of the Icnsing galaxy. Recent work bv lYoo et al. (2005|) 
suggest that the position of the lensing galaxy is an ap- 
propriate approximation of the center of the lensing po- 
tential. For lens systems with more than a single source 
position (i.e., jet sources), quadruply imaged systems will 
always be overconstrained without reducing the number 
of parameters. Note that in all sections, we restrict our- 
selves to four image lenses. 

3. INTRINSIC ASTROMETRIC PERTURBATIONS 

In this section we investigate the intrinsic perturba- 
tions generated by substructure: the position difference 
between images generated by a parent halo alone and 
images generated by a parent halo with substructure. 
While observationally we are interested in the residuals 
relative to the best-fit model, comparisons to the intrinsic 
perturbations allow us to determine to what extent per- 
turbations by substructure are degenerate with changes 
in the macromodel parameters of the best-fit lens. 

Our initial expectations for intrinsic perturbations are 
colored by our understanding of single perturbers. Con- 
sider the case of adding a single perturber to a smooth 
macromodel; it is clear that the astrometric perturba- 
tion will scale with the the size and radial position of 
the perturber, so that larger and more centrally located 
subclumps produce larger perturbations. In the case of 
substructure distributions, however, it is unclear if the 
most massive/most centrally located substructures will 
dominate the astrometric perturbations. For instance, 
the steepness of the mass function means that there are 
many more small perturbers than large, so if the for- 
mer act cooperatively they could in principle generate 
a large perturbation. Conversely, since oppositely posi- 
tioned perturbers generate equal and opposite perturba- 
tions, the net effect of a large number of substructures 
may cancel out ensuring that rare, massive substructures 
dominate the position perturbation of the images. 

Figure [3] shows the intrinsic perturbation distribution 
for each of the fiducial substructure models considered; 
each system is quadruply-imaged and the residual for 
each image position is calculated. We note here that 
in this figure - as in all of our subsequent figures in- 
volving histograms of residuals - shows the results in a 
logarithmic scale. Testing the entire subclump distribu- 
tions, we find that the residual distributions all have very 
large peak perturbations, > 10 mas. Interestingly, our 
simulation-derived substructure models result in residu- 
als that are comparable to or larger than those of the 
semi-analytic models, demonstrating that the intrinsic 
astrometric perturbations are not necessarily dominated 
by nearby clumps. In addition, since the simulation mod- 
els have extremely few or no substructures near the Ein- 
stein radius of the lens and therefore the perturbers must 
be located further away, we can infer that position per- 
turbations of different images in any lens configuration 
may be strongly correlated. A final interesting quality of 
our residual is that the two simulations with the largest 
peak residuals are also the two models that have very 
massive substructures. This suggests that rare, massive 
clumps may cause larger perturbations than the more 
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Fig. 3. — Histograms of the residuals between observed images 
in milliarcseconds for systems with substructure and without sub- 
structure for quadruply-imaged single point source systems. The 
panels show the z = 0.5 simulation distribution (top-left) dynami- 
cally young semi-analytic substructure distribution (top-right) , the 
dynamically old semi-analytic distribution (bottom-right), and the 
2 = simulation distribution (bottom- left). 

abundant smaller clumps. However, even in the mod- 
els where no such massive substructures are present, the 
astrometric perturbations of the images is still consider- 
able. 

4. MODELED RESIDUALS 

In general, the intrinsic perturbations just discussed 
are at least partly degenerate with macromodel param- 
eters. For instance, substructure can change the macro- 
model by adding to the total projected mass distribution 
and increasing the Einstein radius, 6; by changing the el- 
lipticity or orientation of the macromodel; and by adding 
external shear to the potential. For example, the pertur- 
bations from a single subclump placed far from the lens 
are degenerate with an external shear. We expect then 
that the modeled residuals, residuals between the ob- 
served image positions and the image positions predicted 
by the best fit smooth macro-model to the lens, will be 
smaller than the intrinsic residuals discussed above. 

We test the extent to which the intrinsic perturbations 
can be accommodated with the macromodel and whether 
the remaining perturbations (i.e., the perturbations rel- 
ative to the best-fit model) are large enough to be de- 
tectable. We consider only quadruply-imaged systems 
and test two possible lens modeling scenarios: (1) a sin- 
gle point source system modeled with the 7 parameters 



mentioned earlier (&, q, 9q,"f, 9. 



, with the center of 
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the potential fixed) and (2) a jet source, approximated by 
two source positions separated by 10 mas and modeled 
by 11 parameters, including the center of the potential 
and two source positions. In jet source cases, we refer 
to each pair of images produced by the source positions 
together as an image of the jet, and refer to the image of 
each individual subsource as a subimage. 

4.1. Automated Lens Fitting 

Before we move on to results, a discussion of one of 
the difficulties in the automated lens fitting procedure 
is required. The lens fitting is done through a simple 
downhill simplex algorithm. If our initial guess for the 
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Fig. 4. — Histograms of the residuals between modeled and 
observed images in milliarcseconds for single point source systems 
using 7 fitted parameters and fixing the center of the lens potential 
to (x,y)=(0,0). Only four image systems are plotted. In each 
panel the solid line indicates images with Gaussian errors and no 
substructure, while the dotted line is the same, where the initial 
guess for the Einstein radius is underestimated by 10%. Clockwise 
from the top-left panel, Gaussian errors with a standard deviation 
of 0.1, 1, 10, and 100 mas are shown. We can see that a tail 
of large residuals is an artifact of our fitting procedure, though 
the peak perturbation scale is robustly determined. Consequently, 
in all future plots of residual distributions we focus only on the 
scale at which the distribution peaks. The tail of high residuals 
comprises about 10% of the total number of sources considered for 
each model (si 10000). 

best-fit model is not in the same valley as the true 
best-fit model, our resulting formal best-fit model can 
be quite far from the true best-fit smooth macromodel. 
We illustrate this point in Figure ID For this figure, we 
generated lenses using the smooth host halo mass distri- 
bution only (no substructure) , and then added Gaussian 
perturbations of 0.1, 1, 10, and 100 mas to the image 
positions in successive panels. The resulting lenses were 
then fit starting from two different initial guesses for the 
best fit model. For the first, shown in Figure [Has a solid 
line, we used the original mass distribution used to gen- 
erate the lens in the first place (i.e., the "correct" smooth 
mass model). For the second, shown with the dotted line 
in Figure m we simply lowered our initial best guess for 
the Einstein radius by 10%. In all cases, only four image 
systems are used and the residual from each image is cal- 
culated. As is evident from Figure 21 in the first case the 
lens is always fit correctly, and the residual distributions 
are sharply cut around the scale of the deviation, exactly 
as we would expect. On the other hand, for the second 
case we find that there is a large, unphysical tail of high 
residuals. This tail was found in all of our substructure 
runs, and accounts for about 10% of the systems fit. 
This relation is true for Gaussian perturbations regard- 
less of the size of the perturbations. 

In the remaining sections of our work, we ignore this 
tail as an artifact of our lens fitting algorithm. More 
importantly, the peak of the residual distribution is cor- 
rectly recovered by our algorithm, and the location of the 
peak in our substructure runs has physical significance. 

4.2. Substructure Degeneracies with Macromodels 

We present the residual distribution for the single 
source and jet source scenarios using realistic substruc- 
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Fig. 5. — Histograms of the residuals between modeled and 
observed images in milliarcseconds for single point source systems 
using 7 fitted parameters and fixing the center of the lens potential 
to (x,y)=(0,0) (solid) and for jet source systems using 11 fitted 
parameters (dotted). Only four image systems are plotted. The 
panels show the z = 0.5 simulation distribution (top-left) dynami- 
cally young semi-analytic substructure distribution {top-right) , the 
dynamically old semi-analytic distribution (bottom-right), and the 
z = simulation distribution (bottom-left) . 

ture distributions in Fig. [Sj Here, the modeled residuals 
are small and, for the most part, significantly smaller 
than the intrinsic residuals discussed in the previous sec- 
tion, with the peak perturbations occurring between .1 
mas and 10 mas. The decrease in the size of the peak 
perturbations with respect to the intrinsic perturbations 
reflects the fact that there is significant degeneracy be- 
tween the subclumps and the macromodel. 

Comparing the different substructure distributions, we 
see that for the single source case, the semi- analytic sub- 
structure models produce larger perturbations than the 
simulation substructure models. This difference is evi- 
dence for the importance of having substructures within 
the neighborhood of the Einstein radius of the lens in 
order for the system to exhibit observable astrometric 
perturbations. It is also worth noting that the z = 
simulation substructure model had the largest intrinsic 
position residuals and has the largest fraction of poor 
modeled fits - demonstrating that a massive but distant 
clump changes the macromodel enough to be more dif- 
ficult to find a best-fit model. The acceptable modeled 
residuals, however, are comparable to that of the z — 0.5 
simulations, showing that those perturbations induced 
by a massive but distant clump are degenerate with an 
external shear. 

In addition to substructure model differences, there are 
significant differences between single point source sys- 
tems, where we hold the center of the lens potential fixed, 
and jet source systems, in which the center is allowed to 
float. When jet sources are used, modeled residuals are 
smaller than the corresponding residuals in the single 
source case. This result is contrary to our naive expecta- 
tions. One could imagine taking each set of four subim- 
ages, and fitting each individually with a mass model. 
Since the best-fit models for each set of four subimages 
will in general differ, a single mass model for all eight 
images should result in larger residuals. The resolution 
to this problem becomes apparent when we realize that 
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Fig. 6. — Histograms of the residuals between modeled and 
observed images in milliarcseconds for images with Gaussian errors 
with a standard deviation of 1 mas and no substructure. For jet 
sources, each subimage is perturbed separately in the dot-dashed 
line but both subimages are perturbed by the same amount in 
solid line. In both panels, the shaded, dotted histogram shows 
the distribution when employing single point source systems, with 
the center of the lensing potential fixed. Top: The center of the 
lensing potential is held fixed for the jet source systems. Bottom: 
The center of the lensing potential is fit for the jet source systems. 



ber of a subimage pair is perturbed in the same manner, 
however, the distribution of residuals changes dramati- 
cally, and in particular the peak residual is an order of 
magnitude lower than what we found when the center of 
the potential was held fixed. Thus, our observed resid- 
ual distribution from the substructure lenses indicates 
that the position perturbations generated by substruc- 
tures are coherent on scales at least as large as the as- 
sumed extent of our jet source. 
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in fitting jet sources we have allowed the center of the 
potential to float, so the added freedom could result in 
a better fit. Now, in practice, we find that this added 
freedom does result in better fits, but this didn't have to 
be the case, i.e., substructure position perturbations did 
not have to be degenerate with a change in the position 
of the lens's center a priori. 

We illustrate this point in Figure [6l The figure is pro- 
duced by adding 1 mas Gaussian errors to images for 
both single point source systems and jet source systems 
as lensed in the absence of substructures (i.e., we use 
only the smooth halo component to generate the artifi- 
cial lenses). In addition, we test two methods of adding 
the Gaussian errors to the jet source systems; in one, 
each subimage is perturbed independently of the other, 
while, in the other, each image is perturbed separately, 
but the two subimages of the each lensed image are per- 
turbed in the same way. We then proceed to fit the 
artificial lenses. As before, single point source systems 
arc modeled by holding the center of the potential fixed. 
Jet source systems, on the other hand, are modeled both 
with the center of the potential fixed and with the center 
of the potential allowed to float. For the former case (top 
panel), we find that, just as we would expect, the single 
source case (filled, dotted histogram) results in smaller 
residuals than the jet source case (dot-dashed line) when 
each subimage is perturbed independently. When each 
subimage pair is perturbed in the same way, the result- 
ing histogram (solid line) is essentially identical to that 
of the single point source case, again in agreement with 
our expectations. Turning now to the bottom panel of 
Figure [6] in which the center of the potential is allowed 
to float, we find that this added freedom leaves the resid- 
ual distribution of the jet source unchanged when each 
subimage is perturbed independently. When each mem- 



FlG. 7. — Histograms of the residuals between modeled and 
observed images in milliarcseconds for single point source systems 
using 7 fitted parameters and Gaussian priors on the center of the 
lens potential of 5 mas (dotted) and 1 mas (solid). Only four image 
systems are plotted. The panels show the z = 0.5 simulation dis- 
tribution (top-left) dynamically young semi-analytic substructure 
distribution (top-right), the dynamically old semi-analytic distribu- 
tion (bottom-right) , and the 2 = simulation distribution (bottom- 
left). 

Given the importance of fixing the center of the poten- 
tial, an interesting question to ask is how well must the 
center of the lensing potential be known for astrometric 
perturbations to be sizeable. We address this question 
in Figure [7l where we test the residual distributions of 
lenses with substructure for single point source systems 
where the center of the potential is constrained with a 
Gaussian prior. As can be seen from the figure, a Gaus- 
sian prior of 5 mas does not result in much improvement 
relative to the case with no prior (the jet source case), 
while a prior of 1 mas results in a residual distribution 
similar to that obtained when we fixed the center of the 
potential (our fiducial single point source case). Thus, we 
expect that the center of the potential must be known to 
w I mas accuracy for substructure perturbations to be 
non-degenerate with an allowed shift of the center of the 
lensing po t ential . As discussed previously, the results of 
lYoo et all ()2005D suggest that the assumption that the 
lensing galaxy represents the center of the lens poten- 
tial holds to within 5 to 10 mas. Detailed model fits 
of Hubble Space Telescope (HST) data shows that lens 
galaxy astrometry down to 2 mas are achievable (e.g., 
Ilmpev et aLlllQQSli iLehar et al.ll2000[ ). 

4.3. Massive Substructures Near the Einstein Radius 
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Fig. 8. — Histograms of the residuals between modeled and ob- 
served images in milliarcseconds for our alternate projection of the 
z = simulation substructure model that projects a large sub- 
clump near the center of the system. Only four image systems are 
plotted. Top: Single point source systems using 7 fitted parameters 
and fixing the center of the lens potential are plotted in the solid 
line and jet source systems using 11 fitted parameters are shown 
in the dotted line. Bottom: Histograms of the residuals between 
modeled and observed images in milliarcseconds for single point 
source systems using 7 fitted parameters and Gaussian priors on 
the center of the lens potential of 5 mas (dotted) and 1 mas (solid) . 
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Fig. 9. — Histograms of the residuals between modeled and 
observed images in milliarcseconds for the dynamically old semi- 
analytic substructure model. Only four image systems are plotted. 
Single point source systems using 7 fitted parameters and fixing 
the center of the lens potential are plotted. The fiducial result is 
shown by a solid line. Increasing the concentration by 2 (at fixed 
mass) is shown by a dotted line and increasing the concentration 
by 10 is shown by a dashed line. 

One additional possibility we consider is what happens 
when an extremely massive substructure projects near 
the Einstein radius of the lens. This question is rele- 
vant not only because one expect such cases to result in 
the largest residuals, but also because, observationally, a 
large fraction of the current lens sample is seen to have 
luminous satellites projected near the Einstein radius. 
While luminous satellites can be directly included in the 



lens model, it is still interesting to ask what residuals 
would a similar dark substructure produce. 

To address this question, we have chosen an an alter- 
nate projection of the z = simulation model in which 
the most massive halo substructure gets projected to 
within ~ 2" of the halo center. We choose the z = 
simulation model because this model contains the largest 
subclump among our four substructure realizations with 
a mass of about ^ 7% of the host halo mass. A dark sub- 
halo of this size is unlikely, but it is interesting to test the 
extremes necessary for large astrometric perturbations. 

The resulting residual distributions for this substruc- 
ture model are shown in Figure [5] for both a single point 
source (top panel, solid line) and a jet source (top panel, 
dotted line). The same basic qualitative features as seen 
in the fiducial cases are found in this projection. The 
peak residual for the single source case is ~ 10 mas and 
drops to « 1 mas for the jet source case. The total 
amount of substructure within 2" of the lens center in 
this projection is significantly larger than those found 
for the fiducial substructure models, but the peak resid- 
uals are only somewhat larger. We also show in the bot- 
tom panel the residual distributions obtained using 1 mas 
(solid) and 5 mas (dotted) Gaussian priors on the cen- 
ter of the lensing potential and using only a single point 
source. We again find that the center of the lensing po- 
tential must be known to within about 1 mas for the 
position residuals to be sizeable. 

4.4. Dependence on the Concentration of Suhhalos 

In the previous sections, we have preferentially chosen 
substructure realizations that would maximize the num- 
bers of subhalos projected near the center of the lens - us- 
ing semi-analytic substructure models and choosing the 
line-of-sight along the major axis - thereby maximizing 
the estimated signal from substructure. For distributions 
with significant numbers of subhalos projected near the 
Einstein radius, the peak of the residual distribution is 
a few milliarcseconds. As discussed previously, however, 
the subhalos in both the simulation-based substructure 
models and the semi-analytic substructure models could 
be biased to lower concentrations. Higher concentrations 
may lead to larger deflections but may require better 
alignment with images. In Figure [51 the residual dis- 
tribution for the fiducial dynamically old semi-analytic 
substructure realization is compared to the same realiza- 
tion where the concentration of each subhalo is doubled 
while the mass and tidal radius of each subhalo is held 
fixed and to the realization where the concentration is 
increased by an order of magnitude. Here, it can be seen 
that the peak of the residual distribution is lowered when 
the concentration of clumps is increased, so, in fact, our 
estimates likely represent the largest values for astromet- 
ric perturbations. 

4.5. Halo-to-Halo Variation 

The fiducial substucture distributions show that the 
size of modeled residuals increases with a greater amount 
of subhalos projected near to image positions. Halo-to- 
halo variation, however, could swamp this effect. A full 
statistical analysis of halo-to- halo variation would require 
a more extensive understanding of both observational bi- 
ases and substructure model parameters. In this section, 
we test the possible range of the scatter using 48 different 
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Fig. 10. — The peak of the residual distributions for 4-image 
single point source lens systems for 48 semi-analytic substructure 
models using a random hne-of-sight. The average residual peak 
size is shown as square points, enclosed within the minimum and 
maximum peak residuals values. Top: The average peak residual 
compared to the number of subhalos within 2" of the halo center. 
Center: The average peak residual compared to the total mass 
in subhalos within 2" of the halo center. Bottom: The average 
peak residual compared to the mass of the largest subhalo in each 
distribution. 

substructure distributions drawn from the semi-analytic 
models discussed previously and including the two fidu- 
cial semi-analytic models. Testing single point source 
systems with a random line-of-sight, these semi-analytic 
substucture models have between zero and fifteen sub- 
halos within 2 arcsecs of the center of the halo. Most of 
the substructure realizations contain between three and 
five clumps within 2 arcsec of the center of the halo. The 
number of models with a particular number of nearby 
clumps cannot be ascribed any particular significance, 
however, a relationship between the number of nearby 
clumps and the mean of the peak of the residual distri- 
bution may be seen in the top panel of Fig. (TUl The 
range of residuals peaks in the bins is also shown in this 
figure, where, for the most populated bins, the scatter is 
significant, but not as large as the total range of peak 
distributions. For all the models, the residuals fall be- 
tween 0.1 mas and 30 mas. A similar correlation is found 
between the average peak residual and the total mass of 
subclumps within 2 arcsec. For comparison, we also plot 
the average peak residual compared to the mass of the 
largest subhalo. Here, no clear correlation can be ob- 
served. 

4.6. Discussion 

Our results indicate that systems with massive sub- 
structures within the area of the Einstein radius result 
in the largest residuals, of a few mas, which are an or- 
der of magnitude larger than residuals for substructure 
distribution with no nearby clumps. Just as importantly, 
however, we have found that without milliarcsecond con- 
straints on the center of the lensing potential, astromet- 
ric perturbations from substructures are sufficiently de- 
generate with changes in the position of the lensing po- 
tential so as to reduce modeled residuals by an order 
of magnitude and making the residuals extremely diffi- 



cult t o detect on the b asis of astrometric perturbations 
alone. lYoo et all ()2005f ) have suggested that the assump- 
tion that the lensing galaxy represents the center of the 
lens potential holds to within 5 to 10 mas. Astrometry of 
less than a few milliarcseconds for lensed galaxies may be 
achieved using detailed model fits of Hubbl e Space Tele- 
scope (HST) data (e.g.. Ilmpev et al.l 119981 : iLehar et al.1 
[2000h . 

Interestingly, iBiggs et"an (|2004[ ) have presented a four 
image jet source system, CLASS B0128+437, in which 
reproducing the observations with a smooth mass model 
appears to be difficult, even while allowing the center of 
the lensing potential float. More specifically, the largest 
position residuals they find are of order a few mas. Such 
large perturbations are comparable to the highest peak 
residuals we found, and suggest a large substructure mass 
fraction at ra dii near the Einste i n rad ius. The errors 
in the original [Palal fc Kochane"^ (|2002[ ) result are large 
enough to include both small and large substructure mass 
fractions (0.6% to 7% within the 90% confidence inter- 
vals) . In simulations, iMao et al] (|2004D found that the 
substructure mass fraction at radii near the positions is 
smaU, ~0.5% - as small as the amount of substructure 
found in our projected simulation halos which have small 
astrometric residuals. In addition, studies of cusp rela- 
tion anomalies with ray-tracing of simulations of galac- 
tic halos with substructure have implied that there in 
not enough substructure in simulat ions to account for 
the observed level of cusp anomahes (iBradac et al.|[2004l : 
lAmara et al."2006': 'Mac cio et al.ll2006i ) 

Yoo et al. (2005) measured the displacement between 
the observed position of various lensing galaxies and the 
position of the best fit lensing model. Their finding of 
displacements between the two of order 5 to 10 mas im- 
ply that the center of the lens potential and the posi- 
tion of the lensing galaxy are at least roughly coincident, 
but they also suggest that the displacements may in- 
dicate an alternative way to measure the substructure 
mass fraction. Displacements of this size are consistent 
with the typical displacements we observed, modeling 
our clumpy lenses with smooth macromodels. However, 
when we allowed the center of the potential to float using 
jet sources, we found that that even when there were no 
substructures or an extremely small substructure mass 
fraction within 2" of the center of the lens, the center of 
the potential was still displaced by a few to tens of mil- 
liarcseconds. This suggests that such displacements are 
sensitive to substructure far from the Einstein radius of 
the lens, and hence interpretation of such observations 
as limits on the substructure mass fraction in the cen- 
tral regions of lenses might be subject to some caveats. 
Further investigation of whether the substructure mass 
fraction can be robustly estimated from mis-alignment 
of the observed position of the lensing galaxy with the 
best fit location for the center of the potential is clearly 
warranted. 

5. DETECTING SUBSTRUCTURE IN JET SYSTEMS 

Given the difficulties in using lens modeling to probe 
substructure distributions, unambiguous detections of 
substructure in individual lenses might be more robustly 
achieved through local measures of astrometric pertur- 
bations. By offering more observational constraints, jet 
source systems seem particularly promising avenue of re- 
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search. 

5.1. Astrometric Signatures of Substructure 

Given a perfectly straight jet source small enough that 
lensing can be locally linearized, its images will also be 
perfectly straight. If one were to observe a lensed jet 
with two images, one straight and one with a kink, one 
might be tempted to conclude that the kink must have 
been introduced by a local small scale perturbation to the 
lensing potential. However, large kinks may also arise in 
the absence of substructures since small deviations from 
linearity in the jet can be greatly amplified when the lens 
mapping is nearly singular.^ In general, then, without a 
priori knowledge of the linearity of a jet source, some 
degree of modeling is required to detect substructure. 

5.2. Local Lens Modeling 

In this subsection we consider one possible method for 
detecting astrometric perturbations using purely local in- 
formation from the images (i.e., eschewing global lens 
modeling). We consider lensed images of sources where 
multiple source subcomponents are resolved and where 
each of the subcomponents is multiply imaged. More- 
over, we assume the images are in a fold configuration, 
and focus our attention on the close pair of images. As- 
suming the size of the source is smaller than any length 
scale associated with the lensing potential, the mapping 
between the source and each of its images may be lin- 
earized. Letting and A- denote the inverse magnifi- 
cation tensors describing the mapping from the positive 
and negative parity images, respectively, onto the source, 
the two images must themselves be related via a linear 
transformation, 6x^^ — U ■ (5x_, where U — {A+)~^A^ 
and (5x± denotes the image position vectors in any co- 
ordinate system chosen such that the origins in the lens 
plane maps to a single origin point in the source plane. 
Moreover, the assumption that the images are in a fold 
configuration puts a constraint on the form of the ma- 
trix U , so deviations on the form of the matrix U from 
its expected structure may signal the presence of sub- 
structures on scales comparable to or smaller than the 
separation between images. The goal of this section is to 
investigate this possibility quantitatively. 

Operationally, we proceed as follows. Given a jet with 
three subimages, we define the best-fit linear transfor- 
mation U by minimizing the total residual (Ar) defined 

by 



(Ar) = 




)~U-ix 



(5) 

where the sum is over the position vectors, x^^' is the 

image position of vector z, Xj^-* is the image position cho- 
sen as the origin of the coordinate system, and N is the 
number of terms in the sum. For three jet subcompo- 
nents, = 2, so fitting is overconstrained if the matrix 

® While bent radio jets do not, in general, imply the presence of 
substructure, some lens configurations can be described that would 
unambiguously signal the presence of localized perturbations. For 
instance, if two jets in a fold configuration — where two of four 
lensed images are close together — are bent in the same sense, one 
would associate the same parity to both images, violating a generic 
prediction of smooth lensing potentials. 



U has three free parameters or less. Here, we constrain 
the form of the relative distortion matrix U using the 
well known fact that for smooth lenses fold images are 
expected to have equal and opposite magnifications as 
the image separation goes to zero. We thus set the con- 
dition det(U) = — 1; violations of this condition may be 
indicative of substructure. It is worth noting that this 
constraint is coordinate independent, so fitting can be 
done in any appropriate, observationally-defined coordi- 
nate system without any loss of generality. For com- 
pleteness, appendix [^explicitly performs a Taylor series 
expansion of the lensing potential around a specific point 
along the lens's critical curve and in a specific coordinate 
system, though we found this more detailed analysis to be 
less sensitive to substructure than the simpler approach 
presented here. 

The properties of the distribution of the residual Ar 
also depend on the assumed properties of the jet source. 
The magnitude of the residual scales with distance be- 
tween jet subcomponents, which we take to be 1 mas 
in our examples. More interestingly, we will find that 
the distribution of Ar depends on the jet bending angle, 
the angle between the two relative position vectors of 
the source subcomponents. Thus, a zero degree bending 
angle corresponds to a perfectly straight jet. 

Figure [TT] shows the (Ar) distribution while enforcing 
det{U) — —1 for various jet bending angles. Solid lines 
are obtained in the no substructure case, and the cross- 
hatched histograms correspond to the case of the dynam- 
ically old semi-analytic substructure model. We can see 
that for straight or nearly straight jets, the substructure 
and no substructure cases are indistinguishable. As the 
jet becomes non-linear (the bending angle increases), the 
Ar distributions for the substructure case starts to ex- 
tend to higher Ar values, implying that nonlinear jets 
may in principle detect substructure through this local 
test, but only if the substructure density is high enough. 
Indeed, the simulation substructure models are always 
indistinguishable from the no substructure case, which 
is not unexpected given the lack of substructures in the 
central regions of the halo. 

Our result in Figure [Tl] might seem surprising: linear 
jets are ineffective at detecting substructure without the 
full lens modeling. As suggested strongly by section 
the simple picture of substructure producing a bend in 
a single image while leaving the other untouched seems 
overly naive. Moreover, jets that are bent to begin with 
can detect substructure. We interpret the bending angle 
dependence as follows: for a perfectly linear jet, obser- 
vations of the jet constrain the lensing distortion along 
the single jet axis only, so the correct amount of distor- 
tion can be reproduced in the substructure realizations 
regardless of the det{U) — —I constraint. For bent jets, 
however, the distortion along two different axis is probed, 
so both eigenvalues of the relative distortion matrix can 
be constrained by observations. Since the det(J7) = — 1 
condition specifies a relation between these two eigenval- 
ues, it is not surprising that only when both eigenvalues 
are resolved can the effects of substructures be discerned. 

6. LUMINOUS SATELLITES 

In the previous sections, wc have focused on cases 
where the substructure producing astrometric perturba- 
tions is dark. However, luminous satellites can perturb 
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Fig. 11. — Distribution of (Ar) for different jet bending angles. 
Clockwise from top-left, the bending angles are 0, 30, 60, and 90 
degrees. The thick solid line shows the distribution for the smooth, 
no substructure model, while the cross-hatched region shows the 
distribution for the old, semi-analytic substructure model in the 
0, 30, 60 degree panels and increasing in deviation from the no 
substructure case In the 90 degree panel, an additional shaded, 
dotted region shows the distribution for the z = simulation sub- 
structure model, which does not appreciably deviate from the no 
substructure case. 

image positions as well, and in a number of strong lens- 
ing systems, bright satellites have been detected near 
the lensed images. Such systems are important for sev- 
eral reasons. First, and most obviously, much more can 
be learned about the subhalo if its position is known: 
since the astrometric perturbation Sa (x m/r for a satel- 
lite of mass m and distance r, a given perturbation 
can be produced by a nearby low-mass subhalo or a 
distant high-mass subhalo. Additionally, by comparing 
the lensing properties of luminous satellites to those of 
dark subhalos, we can address why are some subhalos 
are dark and some are luminous. Gravitational lens- 
ing can help resolve this question by directly measuring 
the masses associated with satellites falling near lensed 
images. We illustr ate this by co nsiderin g two systems , 
MG J0414-I-0534 ( Trotter et all 2000: r'os et all [2f)M) 
and MG J2016+112 ( Koopmans et al.ll200l ). which both 
have lenses at redshift 2; ~ 1. A key advantage of 
these systems over other lenses with known satellites 
(e.g. RXJ0911-I-0551) is that both have been observed 
with high-resolution VLBI imaging, giving astrometry 
at the sub-milliarcsecond level. We infer the satellite 
masses in these systems by modeling the lensed image 
positions only, ignoring the image fluxes. We hold fixed 
the position of the main lens galaxy and the satellite at 
the observed locations relative to the images, and model 
the mass distribution using a singular isothermal ellip- 
soid for the main lens and singular isothermal sphere for 
the satellite, also allowing for arbitrary external shear. 
We run a grid of models for many values of the velocity 
dispersion of the satellite, optimizing over all other pa- 
rameters, and define the likelihood as L = exp(— x^/2). 
For MG J0414-I-0534, astrometric d a ta and uncertain- 
ties were taken from iTrotter et al.l (l200Gf) . while for 
MG J 2016-H12 these were taken from iKoopmans et al.l 
(|2002f) . The satellite in MG J0414-t-0534, referred to as 



Object X, has apparent magnitude R = 24.6 correspond- 
ing to absolute magnitude of -19.4 for a ACDM cosmol- 
ogy with = 0.27 and h = 0.71. The satellite Gl 
in MG J2016-I-112 has / = 24.6, corresponding to an 
absolute magnitude of -19.6. 

The results of this procedure are plotted in Figure [121 
For MG J0414-I-0534, the velocity dispersion a of the 
satellite Object X was found to be 81 < cr < 102 km/s 
at 95% confidence. For MJ J2016+112, satellite Gl was 
found to have velocity dispersion 87 < a < 101 km/s 
at 95% confidence. Note that these correspond to quite 
large masses. For comparison, the Milky Way Galaxy's 
largest satellite, the LMC, has an equivalent velocity 
dispersion a « Vrot/V^ ~ 50km s^"'^ and a y-band in- 
tegrated magnitude of ~ —18. While the satellites in 
these systems are apparently rather massive, note that 
the lensing constraints were of sufficient quality to detect 
masses nearly an order of magnitude smaller. 

0.14 1 .... I .... I .... I ... . 
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Fig. 12. — Lensing constraints on satellite masses. The solid and 
dashed curves, respectively, show likelihood as a function of veloc- 
ity dispersion a for the satellites in lens systems MG J2016-1-112 
and MG J0414-I-0534, where the satellite mass distribution is mod- 
eled as a singular isothermal sphere. 

7. SUMMARY AND CONCLUSIONS 

Astrometric perturbations from substructure lensing 
offers a new avenue for probing distributions of DM sub- 
structure in galactic halos. Previous work in this field has 
mostly pr obed the effects of single perturbers on image 
positio ns (jMetealf & Madaulj200l|; llnoue fc Chiball2003l . 
I2005airbh . though Chiba (200% investigated the lens sys- 
tem B1422-I-231 using a simple model in which CDM 
substructures were modeled as point masses - finding 
that typical astrometric perturbations are of order 10 
to 20 mas. Since, observations of the strong lens system 
B0128-I-437 with submilliarcsecond resolution have found 
astrometric residuals to the best fit models to be slightly 
smaller at roughly 5 to 10 mas. Such residuals, however, 
are large en ough to be signifi cant relative to measure- 
ment errors (iBiggs et al.ll2004D . In hght of these finding, 
the question of how large astrometric perturbations by 
CDM substructures are is a very timely one. 

In this paper, we use realistic models of substructure in 
order to create populations of artificial lens and estimate 
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the residuals between modeled and observed image posi- 
tions. Substructure distributions were drawn from DM 
simulations and semi-analytic models in order to bracket 
the expected abundance of dark matter substructures in 
the inner regions of galaxy size halos. Our results can be 
summarized as follows: 

1. In general, astrometric perturbations from sub- 
structure are partly degenerate with the smooth 
macromodel. Intrinsic image position perturba- 
tions from substructure are typically of order > 10 
mas, even when few to no substructures are pro- 
jected near the center of the lens. Modeled resid- 
uals are significantly smaller by about an order of 
magnitude. 

2. Substructure distributions in which clumps project 
near the Einstein radius of the lens appear to be the 
least degenerate with the macromodel parameters 
and, therefore, cause the largest residuals between 
the observed and modeled image positions. When 
the center of the lens potential is fixed, substruc- 
ture models with significant amounts substructure 
projected near the Einstein radius have modeled 
residuals of order a few milliarcseconds. 

3. Substructure astrometric perturbations appear to 
be coherent on scales comparable to the extent of 
typical radio jets. This coherence significantly en- 
hances the degree to which substructure perturba- 
tions are degenerate with a change in the center of 
the lensing potential. Allowing the center of the po- 
tential to be fit as a parameter results in modeled 
residuals that are an order of magnitude smaller 
than would be the case if submilliarcsecond con- 
straints of the center of the lensing potential were 
obtained. 

4. Tests of additional substructure realizations sug- 
gest that ~10 mas may represent the largest resid- 
uals for any realistic substructure distributions 
where the center of the lens potential is fixed. 
These substructure realizations include nearby, 
massive dark subhalos and more concentrated ha- 
los. 

5. Local lens modeling of jet sources can in principle 
detect substructure perturbations provided there 
exist substructures within the Einstein radius of 
the lens in consideration and the jet sources are 
intrinsically bent. 
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6. Using gravitational lensing, it is possible to place 
constraints on the masses of luminous substruc- 
tures of observed gravitational lenses through de- 
tailed lens modeling. These constraints may be 
used to test proposed mechanisms for star forma- 
tion in DM subhalos. 



This work suggests that lens modeling may provide 
an avenue for detecting substructure through astro- 
metric perturbations. Substructure distributions which 
project subclumps near the Einstein radius of the lens 
halo cause the greatest perturbations; however, sub- 
structure catalogs from simulations and semi-analytic 
models suggest that the scale of the perturbations will 
be small, on the order of a few milliarcsecond. How- 
ever, unlike anomalous flux ratios between lensed im- 
ages which h ave been thoroughly inve s tigated in previous 
work s (e.g., |M ao fc Schneidei^ 'IPOS'; 'Metcalf & Madau! 
20011: iDalal fc K ochanck 2002; Kochanck & Dalai 200|; 
Rozo et al.ll2006) . a complete accounting of the effects of 
astrometric perturbations and their observational impli- 
cations has not yet been done; for example, observational 
biases remain unexplored. The results of future studies 
will offer a better understanding of submilliarcsecond res- 
olution observations of strong lens systems. 
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APPENDIX 

TAYLOR SERIES EXPRESSION FOR THE RELATIVE DISTORTION MATRICES OF FOLDS 
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Fig. 13. — The (Ar) (mas) distribution for a sample of realizations of fold image configurations with a 90 degree jet bending angle. The 
thick, solid line represents the no substructure case. The dotted line indicates the dynamically old semi-analytic substructu re di stribution 
and the dashed line corresponds to the 2 = simulation substructure catalog. Left: The first order expansion for U, Eqn. IA6I is tested. 
Right: The zeroth order expansion, Eqn. lAlOl is tested. 



In this a ppendix we pres e nt a brief derivation of the relative distortion matrix U for fold images. We use the 
notation of [Schneider et al] ()1992[ ) throughout. Let then x and y represent position vectors on the lens and source 
plane respectively, and let t/o be the point along the lens's caustic that is closest to the source position it of interest. 
We choose the point yo ^-s the origin of the source plane, and its image on the lens plane Xq as the origin of the 
lens plane. Moreover, we choose our axis so that the inverse magnification tensor A at Xq is diagonal, with the (1, 1) 
component non-zero. 
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With these assumptions. ISchneider et alJ ()1992l equation 6.16) show that one one possible parametrization yc(A) of 
the lens's caustic around yo is given by 

yc(A) = ( M-£^2 )■ (Al) 



Note that A = corresponds to y = 0, that is the source plane origin. Given the source position of interest u, and 
using the fact that, by definition, its distance to the caustic = \u — yc(A)P is minimized at the origin, we obtain 

the constraint equation 

0= f:^l«-yc(A)p 



(A2) 

A=0 



which simplifies to ui — 0. In other words, the coordinate system we have chosen has the property that the source 
position of interest u falls exactly along the y axis. Consequently, from now on we assume u — (0, u). 

Giv en this source position, its corresponding images a;_|_ and a;_ in the lens plane are given by f see ISchneider et al.l 
(|1992D equation 6.17) 

(^±)i=-^ (A3) 
ad 

Phi 

{x±),^±^- (A4) 

where we have defined the Taylor series coefficients a = ipii^ b — (f>^ii2, c — 0^22' ^^'^ ~ '^222- -"-^ is important to 
emphasize that the above results implies that in the coordinate system we have chosen, the image separation vector 
for fold images is along the y axis. Moreover, since the x axis is simply orthogonal to the y axis, the coordinate system 
chosen for the expansion can be observationally determined up to an overall displacement of the origin. 

Returning to our derivation, in terms of the expansi on coefficients, the in verse magnification tensor A{x) at an 
arbitrary lens plane position x near the fold is given by (ISchneider et al.lll992l equation 6.14) 

A(a.)=f ^"^^IrV (A5) 

^ ' \ bXi + CX2 CX\ + dX2 J ^ ' 

so that the relative distortion matrix between the two images of interest is given hy U — {A^)^^A^ where A± — A{x±). 
To lowest order in u, we find 

^ _ ^ 1 + (ei +62) 

where we have defined the quantities 



-fc(l + e2) -l + (ei-e2) 

2c 
1 



u= _v V' r (A6) 



(A7) 



c^-2bd ,/2 

3c^-2bd 1/2 .An^ 

Note that the inverse of this matrix is obtained by simply reversing the signs of ei and £2- The relative distortion 
matrix in the limit u — > is simply given by 

^-{-k -1 ) • (AlO) 

which has a single free parameter to be determined, and det(C/) = —1 as it should. 

Figure [T3| shows the distribution of the average residual ( Ar) given by equation |5| using the zeroth order and first 
order expansions of the matrix U computed above, using the same fold configurations as in Figure [11] for a jet bending 
angle of 90 degrees. The coordinate system chosen for the expansion is approximated by setting the y axis in the image 
plane to be in the direction of the vector between the fold images. As can be seen from the figure, the substructure 
and no substructure residual distributions are indistinguishable. 

It is intriguing that a careful expansion of the lensing potential resulted in a test that is less sensitive to substructures 
than the more simple minded approach taken in the main body of this work. It is possible that this difference simply 
reflects the fact that the expectation det(t7) = — 1 is generic, and not particular to a specific point about which the 
potential is expanded. Since the form of the matrix U explicitly depends on the point along about which the lensing 
potential is expanded, it is possible that a formal expansion about a different origin from the one considered here 
would result in a test of higher sensitivity to substructures. 



